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ABSTRACT 

We investigate the effect of weak gravitational lensing in the limit of small angular 
scales where projected galaxy clustering is strongly nonlinear. This is the regime likely 
to be probed by future weak lensing surveys. We use well-motivated hierarchical scal- 
ing arguments and the plane-parallel approximation to study multi-point statistical 
properties of the convergence field. These statistics can be used to compute the ver- 
tex amplitudes in tree models of hierarchical clustering; these can be compared with 
similar measurements from galaxy surveys, leading to a powerful probe of galaxy bias. 
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1 INTRODUCTION 

Weak gravitational lensing by large-scale density perturbations is likely to be one of the most useful tools for studying the 
statistical properties of the distribution of gravitating matter and the geometry of the background universe (Blandford et 
al. 1991; Miralda-Escude 1991; Kaiser 1992). Studies of lensing in the weak limit usually revolve around the statistics of the 
convergence field k, which is related to the line-of-sight fluctuations in gravitational potential and thus to the distribution of 
Q gravitating mass; for recent technical reviews see Kaiser (1996), Bernardeau (1999) and Mellier (1999). 
^ I Various methods, both analytical and numerical, have been applied to understand weak lensing. Bernardeau et al. (1997) 
C/2 ■ have applied gravitational perturbation theory to compute statistical properties of k. Numerical computations have mainly 
■ focussed on ray-tracing experiments using N-body simulations. Examples of this approach include: Jain & Seljak (1997); 
> ! Jain, Seljak & White (1998, 1999a, 1999b); Premade, Martel & Matzner (1998); van Waerbeke, Bernardeau & Mellier (1998); 
? ""j ' Bartelmann et al. (1998); and Couchman, Barber & Thomas (1998). This method has the advantage that it can treat the small 
' angular scales where perturbative calculations are no longer applicable. It is, however, difficult to probe a large parameter 
space of possible models using numerical experiments and it is consequently important to investigate analytic methods for 
studying the regime where non-linear gravitational fluctuations influence the lensing of light rays. 

In this paper we use a hierarchical model, or ansatz, for the form of the n-point correlation functions of matter in the non- 
linear regime. Such models have a long history (Groth & Peebles 1977; Davis & Peebles 1977; Fry & Peebles 1978; Bernardeau 
& Schaeffer 1992; Szapudi & Szalay 1993) but have recently been shown to generate simple and accurate predictions for count- 
in-cell statistics (Munshi, Coles & Melott 1999a,b), the mass function of collapsed objects (Valageas & Schaeffer 1999) and 
the bias and higher-order correlations of over-dense cells (Bernardeau & Schaeffer 1992; Munshi et al. 1999a; Bernardeau & 
Schaeffer 1999; Coles et al. 1999). Numerical simulations have also demonstrated that the hierarchical ansatz is very accurate 
in predicting clustering properties (Munshi et al. 1999b). It is natural to use this model to study the statistical properties of 
the convergence field involved in weak lensing calculations. 



2 MATTER CORRELATIONS AND CONVERGENCE STATISTICS 

The statistics of n are very similar to those of the projected density field. The formalism for handling higher-order correlation 
functions of projected galaxy clustering was developed by Toth, HoUosi & Szalay (1989) and later applied by, e.g., Gaztanaga 
(1994) and Frieman & Gaztanaga (1994). The relevance and usefulness of these methods to studies of the convergence has 
been pointed out by Gaztanaga & Bernardeau (1998) and Hui (1999). The results we present here are new because they 
concern a different kind of statistical descriptor, the cumulant correlator. We begin this paper, therefore, by developing the 
basic relationship between cumulants (and cumulant correlators) of the convergence field k with the underlying statistics of 
the matter distribution. 

We adopt the following line element for background geometry: 
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dr^ = -cdt^ + a(t)[dx^ + r'^{x)if^\, (1) 

where a{t) is the expansion factor. The angular diameter distance r(x) = K~^^'^ sm[K~^^'^x\ for positive spatial curvature, 
r(x) = (— ^)~^^^ sinh[(— i<')~^''^x] for negative curvature, and r(x) = X for a flat universe. In terms of Ho and fio, K = 
{Q,o — 1)Hq. In the following we consider a small patch of the sky so we can use the plane-parallel (or small-angle) approximation 
to replace spherical harmonics by Fourier modes. The convergence /t at a direction 7 is given by a projection of the (3D) 
density contrast along the line of sight with a weight function a;(x) : 

«(7)= r dxu>{x)S{r{x),^). (2) 
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If all the sources are at the same redshift, one can write the weight function u){x) = S/Aac'^ HQflmr{x)r{xs — X) /"riXs)-, where 
Xs is the comoving radial distance to the source. (This approximation is not crucial and is easy to modify for a more accurate 
description.) Using a Fourier decomposition of 5 we can write 



k(7) 



= J dx<^{x) j -^^e-xp{ixk\\+reki_)5k, (3) 

where we have used fcy and fcx to denote the components of wave vector k parallel and perpendicular to the line of sight. In 
small angle approximation however one assumes that k± is much larger compared to fey , denotes the angle between the line 
of sight direction 7 and the wave vector k. 



2.1 Cumulants 

Using the definitions we have introduced above we can compute the projected two-point correlation function 



(«(70«(7.)) = ['^X.^^ J (I;. eMiOl) ) (4) 

(Peebles 1980; Kaiser 1992; Kaiser 1998) where we defined 1 = r(x)kx, a scaled wave- vector projected on the surface of the 
sky. It is useful to consider the average of this smoothed over an angle ^0, which is often used to reconstruct the matter power 
spectrum P{k) (Jain, Selzak & White 1998). The convergence field smoothed by an angle ^0 is represented by 

Kail) = J d^-yWeoi-j) J dxi^ixJSirix)^), (5) 
where Wbq is the "top hat" window function. We then get 

i-i)= r^x^ f7^p(jz,)wiiWo), (6) 



^^(x) J (27r)2 \r(x)_ 

where W2{x) = 2Ji{x)/x and Ji{x) is a Bessel function; cf. Bernardeau (1995). Similar calculations for the three- and 
four-point correlation functions can be expressed in terms of integrals of the matter multi-spectrum Bp, defined by 

(5(ki)...d(kp)) = Bp(ki,...kp), (7) 

where the kj are constrained to have ^ = 0. First define 



lN[r{x),0o)]= I -^W,{h0o)...J -^W,{l^.,0o)W,{l^eo)B^l-^y...,^^] . (8) 



^From this we obtain 



l3[r{x),0o] (9) 



4>= r dxi^^U[r{x)M (10) 
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and 

r^ix) 

For the derivation of the corresponding hierarchical parameters, 
we need 

In these results we have assumed that the smoothing angle 9o is very small. 



(11) 



(K^) = £'dx jZ^^l,, lNHx),Oo]. (12) 



2.2 Cumulant Correlators 

Cumulant correlators wore introduced by Szapudi & Szalay (1997) as a normalized two-point moments of the multi-point 
correlation functions. Later the concept of two-point cumulant correlators was generalized to multi-point cumulant correlators 
(Munshi, Melott & Coles 1999). The two-point cumulant correlators for galaxies have already been measured in projected 
surveys such as APM by Szapudi & Szalay (1997) and they have been shown to be an important diagnostic of hierarchical 
galaxy clustering models. 

To calculate the two-point cumulant correlators we adopt the additional assumption that the separation between different 
angles are small. This is the regime in which we expect the hierarchical model presented in the next section to be most accurate. 
The general result is: 

(C(7i)<(72)) = £' X+^}^\x) '^X /^■■■J ^7gyF^^^(«i^°) ■ • • W2{lr.+meo) exp[i(Zr + . . . + 1^)9] 

where again = 0, as it is in the following expressions. Relevant low-order examples of this result are: 

(«?(7i)K.(72)) = £' J-§^J -^W,{heo)W,{heo)W,{heo)eMih0^^)B3 ^) ; (14) 

(«?(7i)«.(72)) = P^dx J-^J-^I -^W2iheo)W2il29o)W2{l3eo)eMilsOi^) 

Jl h h h 

r{x)' r{x)' r{x)' r{x) 



(Ks(7l)«s(72)> 



(15) 

r ^"^^ J ^w{heo)w,{i,do)w,{iseo) eMiih + h)e,,) (i6) 
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rix)' r{x)' r(x)' r(x) 



The three-point cumulant correlators can be expressed in the following form 

(C(7l)«?(72)«^(73)> = ^Z+^^r\) '^X J ^■■■J ^^^^W,ih0o) . . . W,{U+r.+reo) (17) 

ll Im + nH 



exp[i(/l -I- . . . -h /m)^12 + i{lm + . . . + lm+n)6l3]B, 



m+n+r 



r{x)'"'' r{x) 



3 TREE MODEL FOR CORRELATION FUNCTIONS 

In deriving the expressions (9)-(17) we have not used any specific form for the matter correlation hierarchy. To go further we 
have to choose a particular form of hierarchical scaling. The most general "tree" model can be written as 

(JV-l) 

Cjv(ri,...rN) = '3iv,c< n ^('■'''■•i)- (^8) 

a, N— trees labellings edges 

A similar hierarchy is developed in the quasi-linear regime in the limit of vanishing variance (Bcrnardeau 1992), but the 
hierarchical amplitudes Qn.u are shape-dependent functions in such a case. In the highly nonUnear regime of interest here, 
there are some indications that these functions become independent of shape (e.g. Scocciamarro et al. 1998). In Fourier space, 
such a result means that the multi-spectra can be written 

B2(ki,k2,k3) = Q[P(ki)P(k2) + P(k2)P(k3) + P(k3)P(ki)] 
B3(ki,k2,k3,k4) = i?a[P(ki)P(ki-Fk2)P(ki-Fk2-^k3)-Fcyc.perm.] 
+Rb[P{ki)P{k2)P{k3) + cyc.perm.] 
S4(ki, k2, ka, k4, ks) = 5a[P(ki)P(ki -|- k2)P(ki +1^2 + k3)P(ki -h ka + kg -h k^) + cyc.perm.] 

-H5'6[P(ki)P(k2)P(ki + k2 + k3)P(ki +k2 + k3+ k4) + cyc.perm.] 

+5c[P(ki)P(k2)P(k3)P(k4) + cyc.perm.] (19) 

(remembering that ^ ki = in all these expressions.) 

Particular hierarchical models differ in the way they fix the various amplitudes for the different topologies (Q, Ra, Rb, 
etc). For example, Bernardeau & Schaeffer (1992) considered the case where amplitudes can in general are factorizable (see 



Munshi, Melott & Coles for a detailed description); see Szapudi & Szalay (1993) for an alternative prescription. We will avoid 
any specific prescription because our aim is to show that the statistics of k can allow one to fix the hierarchical amplitudes 
empirically, as well as constraining the power spectrum and background geometry. We shall ignore smoothing corrections, for 
reasons discussed by Boschan, Szapudi & Szalay (1994). We first define 



(27r)2 yrix) 
r-2(*-i)(x) 



Ct[F{x)] = /' z£^nx)dx. (20) 
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The expression for Ct can be computed numerically for specific models of background cosmology. F{x) denotes the various 
products of K0g and Kg^^^ which appear in the following expressions containing x dependence. The values of hierarchical 
amplitudes (which in general are insensitive to the background cosmology) can be computed from numerical simulations or 
with analytic methods such as hyper-extended perturbation theory (Scoccimarro et al. 1998; Scoccimarro & Prieman 1998). 
The one-point cumulants can then be expressed as 

{«^(7)) = 3Q3C3[k^J 

{/t*(7)> = {12Ra + m)C4[4^] 

{k^(7)> = {60Sa + 60Sb + 5Sc)C5[Ktj (21) 

The first of these expressions was derived by Hui (1999). He showed that this result agrees well with numerical ray-tracing 
experiments by Jain, Scljak & White (1998). The expressions for two-point cumulant correlators are 

(Ks(7i)«s(72)> = aQsCslKSg/veiJ 

(/ts(7l)Ks(72)) = {3Ra+6Rb)C4[KggKei2] 

{k^(-,i)a:?(",2)> = 4_R6C4[KsoKei2] 

^'(7i)'^'(72)) = {12Sa + 6St)C\,[4„^g,2]- (22) 
and the three-point cumulant correlators arc 

{k;s(7i)ks(72)«s(73)> = C5[(45'aKei2K(3i3 -I- 45'6K6ii2Kfli3 + 4S'i,K;ei3«e23)'«flo] 

{/t^(7i)Ks(72)Ks(73)) = C5[[{6Sa + SSb){Kei2i^e23+ 1^613^023) + {SSc + 18Sb+6Sa)Kei2i^ei3]Kla]- (23) 

It is clear from these equations that at any order, measurements of the appropriate cumulant correlators provide enough 
information to determine the hierarchical coefficients of that order and below. Weak lensing statistics can therefore tell us 
directly about the form of hierarchical model appropriate for the underlying density field, just as it can with the distribution 
of galaxies (Munshi, Melott & Coles 1999). 



4 DISCUSSION 

Ongoing optical and radio experiments are intended to provide imaging data for studies of weak lensing. A large smoothing 
angle is required to probe the quasi-linear regime where most existing analytic calculations are valid. To be useful, such 
regions must be at least reach 10 degrees on a side. Existing CCD cameras have typical sizes (0.25° —0.5°). The initial lensing 
surveys will therefore provide information about matter fluctuations in the highly non-linear regime. This is the regime in 
which the model we have presented here is expected to be valid. Ongoing observational programmes in this regime allow the 
direct determination of the hierarchical scaling properties of cosmological density fluctuations. These can be combined with 
similar measurements of galaxy clustering to constrain models of the relationship between galaxies and mass. This is the key 
result of this paper. 

Other studies of small-scale weak lensing have focussed on ray-tracing experiments through N-body simulations. It will 
definitely be very interesting to check our analytical results against such numerical simulations. This will also test different 
approximations wo have used in deriving our analytical results. Results of this kind of analysis will be presented elsewhere 
(Munshi, Jain & White 1999a, b, in preparation). 

Cumulant correlators of orders higher than those we have presented here can be written in a similar way to equations 
(21)-(23). We have not presented them here, however, because sampling variance plays an increasingly important role in these 
quantities at increasing order, and we have not yet developed a method for dealing with this in the context of weak lensing, 
although estimation of finite-volume corrections has been done for galaxies by several authors (Szapudi & Colombi 1997; 
Munshi et al. 1999; Munshi, Coles & Melott 1999b; Hui & Gaztanaga 1998). We also note that the hierarchical ansatz we have 
exploited furnishes a complete prescription of all the higher-order correlation functions of the matter distribution. This allows 
one to estimate the multi-point count probability distribution function for an arbitrary number of correlated volume elements 



(Munshi, Coles & Melott 1999a,b; Bernardeau & SchaefTer 1999). These results can be used to compute the numbers of "hot" 
and "cold" spots in statistics of k in the sky in much the same way as has been done for fluctuations in the cosmic microwave 
background (e.g. Coles & Barrow 1987). The development of our work in these directions will be presented elsewhere (Munshi 
& Coles, 1999, in preparion). 
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